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Abstract 



We have estimated the hmiting branching ratio of the R-parity violating 
(RPV) decay of the hghter top squark, ti — l~^d (/ = e or /i and d is a down type 
quark of any flavor), as a function of top squark mass(mjj for an observable sig- 
nal in the di-lepton plus di-jet channel at the Tevatron RUN-II experiment with 
2 fb~^ luminosity. Our simulations indicate that the lepton number violating na- 
ture of the underlying decay dynamics can be conflrmed via the reconstruction 
of . The above decay is interesting in the context of RPV models of neutrino 
mass where the RPV couplings [X'i^j) driving the above decay are constrained to 
be small ( ~ 10~^ — 10~^). If ti is the next lightest super particle - a theoretically 
well motivated scenario - then the RPV decay can naturally compete with the R- 
parity conserving (RPC) modes which also have suppressed widths. The model 
independent limiting BR can delineate the parameter space in speciflc supersym- 
metric models, where the dominating RPV decay is observable and predict the 
minimum magnitude of the RPV coupling that will be sensitive to Run-II data. 
We have found it to be in the same ballpark value required by models of neutrino 
mass, for a wide range of m^^. A comprehensive future strategy for linking top 
squark decays with models of neutrino mass is sketched. 

PACS numbers: ll.30.Pb, 13.85.-t, 14.60.Pq, 14.80.Ly 



^ spdas <§juphys. ernet. in 
^ adatta@juphys.ernet.in 
guchait@tifr. res. in 



1. Introduction 



The Minimal Supersymmetric Standard Model (MSSM)[1] is a well motivated extension 
of the Standard Model(SM), which is free from several shortcomings of the latter. As of 
now, there is no experimental evidence either in favor of or against it. Unfortunately, the 
mechanism of supersymmetry(SUSY) breaking is not known yet, although several interesting 
suggestions exist [1]. As a result there is no guideline for predicting the mass splitting between 
a SM particle and its superpartner and consequently, there is no information about the 
range of superparticle (sparticle) masses from the theoretical point of view. There are some 
experimental lower bounds from unsuccessful collider searches at LEP [2] and Tevatron Run- 
I [3, 4]. 

Currently the Run-II of the Tevatron (referred to hereafter as Run-II) is in progress. It 
is expected to deliver an integrated luminosity of at least 2 fb^^ per experiment at 2 TeV 
centre of mass energy, which is more than one order of magnitude larger than the acquired 
luminosity in Run-I with center of mass energy 1.8 TeV. However, in view of the existing 
limits on the masses of the strongly interacting sparticles (squarks and gluinos) [2, 3] and the 
rather marginal increase in centre of mass energy, most of the unexplored parameter space 
in this sector is likely to be beyond the kinematic reach of Run-II as well. Yet this is the 
only currently available machine for direct SUSY searches until the LHC starts. 

In view of this, the top squark (the superpartner of the top quark) is somewhat special. 
It may be lighter than the other squarks and gluinos due to several reasons. Firstly, the large 
top Yukawa coupling which controls the evolution of the soft-supersymmetry breaking masses 
of the left and right handed top squarks, f^, in, via the renormalization group(RG) equations, 
tends to reduce these masses [1]. Moreover, because of the large top quark mass, the two 
weak states may mix very strongly leading to a relatively large splitting between the 

two physical mass eigenstates ti,t2 [5] (in our notation m^^ > ^^ifj. Interestingly, the mass 
of the lighter states ii may be even below the top mass. In fact, it is quite conceivable that in 
certain region of SUSY parameter space it happens to be the next lightest supersymmetric 
particle (NLSP), the lightest neutralino x? being the lightest supersymmetric particle (LSP) 
by assumption in most R-parity(-Rp) conserving models. It is, therefore, very important 
to fix up the strategies for isolating the top squark signal for all conceivable decay modes. 
Yet another motivation to look for a light top squark is that it seems to be preferred by 
electroweak baryogenesis [6]. 

In many studies the MSSM is assumed to be a Rp conserving (RPC) theory. The Rp 
is a discrete symmetry imposed on the MSSM to avoid the lepton and baryon number 
violating interactions in the Lagrangian which lead to rapid proton decay. If, however, either 
lepton or baryon number violation is allowed, such catastrophic decays can be avoided. 
This can be achieved by imposing either the so called baryon parity or the lepton parity 
[7, 8] conservation. The resulting theory, called R-parity violating (RPV) MSSM [9], is 
phenomenologically attractive since it has many novel predictions. From the theoretical 
point of view both RPC and RPV versions of the MSSM are on equal footing since both 
require additional discrete symmetries beyond the gauge symmetry. 
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The SUSY signatures arc determined by whether R-parity is assumed to be conserved or 
not. Conservation of R-parity imphes that all SUSY decay chains end up in the Xi: which is 
stable and escapes the detector. Thus a typical SUSY signature is always accompanied by 
some amount of missing transverse energy (^). The R-parity violating interactions [9, 10] 
on the other hand would allow the LSP, which is not necessarily the XiJ to decay into SM 
particles leading to distinct signals [11]. Yet another type of signal is the direct decay of 
sfermions, into two SM particles [10]. In view of the large production cross-section of titi* 
pairs, the class of lepton number violating decays generically denoted by, 

h ^ ltd, (1) 

is an attractive channel for searching RPV interactions at Run-II. Such decays are driven by 
the A^3^ Li Qs Dj term in the superpotential where L, Q and D are respectively the lepton 
doublet, quark doublet and down singlet type superfield and i,j are generation indices. 

In recent times RPV models have attracted special attention as they can provide viable 
models of neutrino mass {nii,). The basic mechanism has been known for a long time[12]. 
The interest in this model was revived after the atmospheric [13] and solar neutrino[14] ex- 
periments confirmed that the neutrinos are not massless. Parameters of these models have 
been constrained by many groups using neutrino data [15, 16]. The actual set of parameters 
( bilinear and trilinear terms [9] in the RPV sector of the superpotential) and their precise 
magnitudes required to explain the neutrino data is model dependent. However some cou- 
plings belonging to the class X'^^j are important ingredients of model building. Considering 
a variety of models it has been shown, e.g., in [16], that the important couplings A^gg, for 
all i, turn out to be generically small ( ~ 10~^ — 10~^, depending on the magnitude of the 
soft breaking parameters in the RPC sector). This is certainly much stronger than the con- 
straints obtained prior to the neutrino data[17]. Thus RPV decays of the top squark driven 
by these couplings may provide an avenue for probing the models of u-mass [18, 19, 20] at 
colliders. 

In this paper we focus our attention on two interrelated topics : 

1. The viability of observing direct top squark decays through the lepton number violating 
channels, Eq.(l), at the upgraded Tevatron collider in a model independent way using 
the event generator PYTHIA [21]. This is an issue important in its own right irrespective 
of models of m^. 

2. The imphcations of observation/non-observation of this decay channel for models of 
neutrino mass. 

The coUider signatures, however, crucially depends on whether the top squark is the 
NLSP or not. If the top squark is not the NLSP and the RPV couplings are as small as 
that required by the models of the neutrino mass, it would dominantly decay via the RPC 
2-body mode with nearly 100% BR, [22, 23] 

h ^ bxt (2) 

where is the lighter chargino, or, if the above mode is not kinematically allowed, via 
the 3-body modes [24], 

tl ^ Mi), biu, bWxl (3) 
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where £ and u are respectively the slepton and the sneutrino assumed to be hghter than ti. 
The decay of the LSP would then be the only signature of RPV interactions. Whether the 
magnitude of the underlying RPV coupling is indeed in the right ballpark or not, can be 
tested in principle, e.g., by measuring the width of the LSP, which may not be an easy task 
at least in the context of Run-II. Of course neutrino masses are generated by specific RPV 
parameters, e.g., by the lepton number violating trilinear couplings A^33, where i is the lepton 
generation index, if the bilinear couplings happen to be small. Thus the decay patterns of 
the LSP may give some circumstantial evidence in favor of/against models of neutrino mass. 
For example, if x? is assumed to be the LSP, then xtxi Xi^Xi production followed by 
decay chains involving the decays 

X? ^ v^hh, Vrhh (4) 

is indicative of an underlying model of neutrino mass [25]. In ref. [25] the prospect of 
observing this signal at Run-II was studied. It was concluded that this signature can be 

probed up to mi/2 = 230 GcV(320 GeV) with an integrated luminosity of 2 fh~^ ( 30 fh^^ 
). Here mi/2 is the common gaugino mass at the GUT scale. It is, however, worth noting 
that since the signal has missing energy carried by the neutrinos, it can be mimicked even 
if R-parity is conserved. For example, the decay X2 ~^ Xi^b, which may have a large BR if 
one of the bottom squark mass eigenstates happens to be lighter at large tan/3, has coUider 
signatures very similar to the decay of Eq.(4). Moreover, the lepton number violating nature 
of the underlying interaction is not obvious since the neutrinos escape the detector. The 
possibility of probing RPV models of neutrino mass through neutralino decays has also been 
considered in [26]. 

The situation is totally different if the lightest neutrahno and the top squark happen to 
be the LSP and the NLSP respectively, a theoretically well motivated scenario for reasons 
discussed above. In this scenario the main RPC decay channels occur via the flavor changing 
neutral current decay mode [22], 

ii ^ cx? (5) 
and via 4-body decay modes with a b-quark, x? and two massless fermions [27], 

ii ^ bx'jf (6) 

( / and /' being a quark-antiquark or pair ) which eventually lead to RPV signals due 
to the LSP decay. Here the key point is that the above channels have widths suppressed 
due to natural reasons and can very well compete with each other [27] or with the RPV 
mode even, if the coupling A-g^ is ~ 10~^ — 10~^. As we shall see in a later section such 
competitions occur naturally over a large region of the MSSM parameter space. In fact if the 
above coupling is much larger then the direct lepton number violating decay mode Eq.(l) 
will occur with 100% BR's. The coexistence of the direct lepton number violating decay 
mode as well as RPC decay modes followed by the LSP decay, is a hallmark of RPV models 
of neutrino mass. Moreover this signal is attractive due to the large production cross section 
of top squark pairs {titi ). In addition, if a signal is observed then the underlying lepton 
number violating interaction can be revealed easily by reconstructing the top squark mass 
as will be illustrated in a later section. The signatures of this decay at the Tevatron was 
first discussed in [28] in the context of high events at HERA [29] . 
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Obviously the presence of competing channels may complicate the search for the top 
squark. For example, if each of the competing modes has BRs substantially smaller than 
100% all of them may be below the observable level in spite of large production rate of titi 
pairs. A complete discussion is not possible without full simulations of all possible signals 
which is beyond the scope of this paper. We shall concentrate on the first task, namely to 
estimate the minimum BR of top squark decay in RPV channel, Eq.(l), required for the 
observation of the signal at Run-II. This will be done in a subsequent section using the MC 
event generator PYTHIA [21]. 

In a recent paper it has been shown [20] that the data from Run-I of the Tevatron already 
restrict the BR of the decay, Eq.(l) to values significantly smaller than 100% in a model 
independent way for a range of top squark mass ( see Fig. 3 of [20]). Assuming specific 
model parameters, which fixes the BRs of all the competing channels, this BR exclusion can 
be translated into upper bounds on the corresponding RPV coupling ( see Fig. 4 of [20]). 
It was found for the first time, albeit for small values of the top squark mass and rather 
limited regions of the MSSM parameter space, that the Run-I data were indeed sensitive to 
magnitudes of these couplings relevant for models of neutrino mass. Since the accumulated 
luminosity of Run-II is more than an order of magnitude larger, we feel encouraged to 
investigate the feasibility of obtaining similar constraints over a much larger region of the 
MSSM parameter space. It may be recalled that in the past Tevatron di-lepton data was 
also used to constrain the squark and gluino masses in the context of RPV SUSY model [30]. 

The possibility of probing the RPV models of neutrino mass via top squark decays was 
also suggested in Ref.[18, 19]. These works, however, differ from ours in several ways. In [18], 
which was the first attempt to confront Run-I data with models of neutrino mass, it was 
claimed that for values of RPV coupling favored by models of neutrino mass, the RPV decay 
of the top squark dominates over the loop decay for m^^ < 150 GeV. In the absence of 
4-body decays this statement is correct for low tan /? only. For high tan j3 the loop decay can 
overwhelm the RPV mode. For low tan/3, on the other hand, the 4-body decay, Eq.(6), not 
considered in [18], may dominate over the RPV decay if A' ^ 10~^ — 10~^. The SU(2) gaugino 
mass M2 also plays an important role in determining the relative strengths of these competing 
modes. All these issues will be addressed in great details in a subsequent section. In Ref. [19] 
the RPV mode, Eq.(l), the loop decay Eq.(5), and several 3-body modes (including those in 
Eq.(3)), were assumed to be the competing channels. However, no detailed simulation was 
carried out to estimate the sensitivity of the data to RPV couphngs. 

We have organized the paper as follows. In Sec. 2, we shall describe our road map for 
obtaining a comprehensive search strategy for top squark and its consequences for models of 
rriiy and briefly review the current status of top squark search, especially when it happens to 
be the NLSP. In Sec. 3, the details of the simulation leading to model independent limiting 
values of the BR (ii — > I'^dj) , where I — e or sensitive to Run-II data will be presented 
as a function of rrif^. In Sec. 4, we use the results of Sec. 3 to obtain upper hmits on RPV 
couplings in specific models and to understand the systematics of the parameter space (i.e. 
delineating the regions where some of the competing modes dominate or several of them may 
coexist). We summarize our results in Sec. 5. 
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2. The road map for linking top squark search with mod- 
els of neutrino mass 



Our first task is to assess the viability of observing the RPV decay, Eq.(l) at Run-II. For 
simphcity we shall as usual assume that the RPV couplings are hierarchical, i.e., one coupling 
of the type X[^j dominates over the others. 

For the sake of definitcness our simulations will be restricted to the mode ti e^dj. This 
decay is triggered by the trilinear RPV coupling X'^^^ L^QsDj term in the superpotential [9], 
where j = 1-3 is a generation index for down type quarks. In order to make our analysis as 
general as possible, we have not employed any particular jet tagging so that the conclusions 
arc approximately valid for any j.^ 

Our conclusions are also approximately valid for the coupling Agg^ Lf^Q-^Dj. A small 
difference may arise due to the difference in the detection efficiencies for e and jj,. However, 
since the leptons are highly central the difference is rather marginal. Our conclusions cannot 
be applied to the signal from the X'^^j L^-Q^Dj term which requires a fresh simulation taking 
into account r detection efficiency We, however, feel that the simplest signal arising from 
the class of decays in Eq.(l) will be sufficiently informative for the first analysis using an 
event generator. 

A systematic search strategy for the top squark or, in the absence of a signal, a compre- 
hensive limit on m^^ in RPV MSSM, therefore depends on several steps. The first step is 
to estimate the model independent minimum value (T{pp ti^i) * (cbr)^ for an observable 
signal as a function of m^^, where ef,r=BR(ti e^dj). Using the well-known formula for 
a{pp ti tl), which is available up to next to leading order (NLO) [31], this bound can be 
translated into a lower limit on observable BR. Rather low values of e^r can be probed for a 
range of rrif^ and rrif^ may be reconstructed with reasonable accuracy at Run-II, as we shall 
see in the next section. 

The observation of the RPV signal alone, though a stupendous achievement in its own 
right, will shed little light on models of neutrino mass(my). As discussed in the introduction 
the simultaneous observation of the signals arising from the RPC decays in Eq.(5) or Eq.(6), 
followed by Xi decay may strongly hint in favor of these models. The observabihty of these 
signals depend on two factors, a) the BRs of the decays involved and b) the acceptance 
efficiency of the cuts in distinguishing the signal from the background. 

Assuming that the dominant RPV decay mode of Xi in RPV models of rrii, is x° —>■ bbui, 
i=l,2,3, the signal resulting from the loop decay (Eq.( 5)) is 

^1 cXi chlvi (7) 

Therefore,tiii* pair production is signaled by jets-|-^ with 4 b-jets. Similarly the 4-body 
decay, Eq.(6) would cascade into 

t^- hxlff bbbujf (8) 

''The most important ingredients of RPV models of are A^33 all of which are constrained to be 
~ 10~^ — 10~^ ( See, e.g, [16] ). However, for dj=h, b-tagging can be efficiently employed to improve 
the signal/background ratio and our conservative conclusions may be further strengthened. ( see Sec. 3 for 
further comments). 
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An excess of £+jcts+^, 2£+jetsH-^ or jets+^ events including several b-jets would, 
therefore, indicate titi* pair production in the framework of RPV SUSY model. The above 
signals are very similar to the ones discussed in [25] although the signal from the decay 
chain in Eq.(8) may have even more b jets. From the results of [25] one has reasons to be 
optimistic that the large number of b jets would provide a visible signal if b-tagging is really 
efficient (a; 50%). 

Full simulations of the above two signals, which is beyond the scope of this paper, would 
lead to the estimated minimum BR of the loop decay and the 4-body decay required for 

observable signals at Run-II as a functions of m^^. These along with the minimum BR for 
observable RPV signal ( estimated in the next section in detail) will provide the basis for a 
model independent approach to top squark search at Run-II in the context of RPV theory 
of nil,. 

If the signal is seen in the RPV channel as well as in one or both of the competing 

channels, one can try to identify the allowed parameter space using the limiting BR and the 
reconstructed m^^. Since the estimates of the limiting BR corresponding to the signals in 
Eq.(7) and Eq.(8) are not available at the moment, a complete job can not be done. However 
the BR of the RPV decays will be discussed in details in Sec. 4 with an aim to understand 
the systematics of the MSSM parameter space vis-a-vis these decays. Outhnes of a future 
comprehensive programme for linking ti decay signals with models of rrii, will be sketched 
with illustrative example in Sec. 4. 

The other important issue is the prospect of unambiguously excluding a range of m^^ if 
no signal is seen. Here one encounters the complications due to possible presence of three 
competing decay modes in a large parameter space. In fact the current mass limits on mf^ 
in both RPC and RPV models are also not free from ambiguities. 

The phenomenology of top squark search Tevatron experiments in different decay chan- 
nels have been studied extensively in both RPV [18, 19, 20] and RPC [32, 33, 34, 35] models. 
The unsuccessful search for the top squark at LEP and Tevatron Run-I experiments in both 
RPC [36] and RPV MSSM [38, 37] have yielded important bounds. Here we shall focus on 
the scenario when top squark is the NLSP. 

It is to be noted that the most stringent limits in RPV as well as RPC models have often 
been derived by employing the model dependent assumption that the top squark decays into 
a particular channel with 100% BR. For example, the most stringent bound in the context 
of RPC MSSM comes from Tevatron Run-I experiments which puts a lower limit on lighter 
top squark mass rrif^ > 119 GeV for m^o =40 GeV. The limit becomes weaker for higher 
value of m^Q, e.g., m^^ > 102 CcV for m^o =50 GeV [36]. In deriving these limits, it was 
assumed that the loop induced, decay Eq.(5) [22], occurs with 100% BR. Apparently this 
assumption is valid in a wide class of models if the ii state happens to be the NLSP. Since 
the production cross section of top squark pairs is dominantly via QCD and depends on its 
mass only, the above limits from Tevatron, therefore, seem to be fairly model independent, 
except for the dependence on m^o. which influences the efficiency of the kinematical cuts. 

However, as has been shown in [27], even if the top squark is the NLSP, its 4-body decay, 
Eq.(6), may indeed compete with the above loop decay or may even overwhelm it in some 
region of parameter space. The above limits, therefore, require revision and new signal via 
the 4-body decay channel should be looked for [35]. 
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The most recent limit on the top squark mass ( m^^ ~ 122 GeV ) in the RPV MSSM 
comes from the CDF collaboration [37] in the decay channel 

h ^r+ + b (9) 

This limit is also derived on the basis of the above model dependent assumption, namely, the 
decay channel in question has a BR of 100%. However, even if the RPV coupling involved 
(A333 ) is assumed to be the most dominant one, the mode may have a BR significantly smaller 
than 100 %. This may happen in various regions of the MSSM parameter space simply due 
to the competition among this decay mode and several RFC modes of top squark, since the 
latter couphngs are invariably present in the theory irrespective of the choice of the RFV 
sector. As discussed in the introduction the competition is of special interest, if the top 
squark is the NLSF and RFV couplings have strengths relevant for the models of neutrino 
mass [15, 16]. In [20] on the other hand the possibility of competition among different decay 
channels were considered. The mass limits obtained in [20] were naturally dependent on 
BR(ti ed )= e^r. For example it was found that rrif^ ~ 200(165) GeV for €6^=1(0.5). 

If no RFV signal is seen at Run-II, any particular m^^ can not be excluded in a model 
independent way. Only the regions of the MSSM parameter space where the BR of at 
least one of the three competing decay modes is above the observable limit will be ruled 
out. On the other hand one can also identify the difficult regions of the MSSM parameter 
space, in the context of Run-II , where all three decay modes have low BRs. The stop 
search at LHC may focus on these regions. In the difficult regions the RFV signals from 
chargino/neutralino production followed by Xi decay [25] appear to be the only possibility 
of probing models of rriu at Run-11 . Thus the top squark decay and the signal of [25] are 
essentially complementary in nature. 

It should also be noted that the above top squark decay signals are important only if the 
top squark happens to be the NLSF, a scenario theoretically very well motivated but not 
inevitable. The signal of [25] on the other hand requires the lighter chargino to be heavier 
than the Xi which is not necessary in RFV models, unless gaugino mass unification [1] is 
assumed. Thus either of the above two signals, Eq.(l) or Eq.(4), or both may be helpful for 
probing RFV signals depending on the MSSM parameter space of interest. 

The limit on the RFV BR in turn can be converted into upper limits on A' in specific 
models with several competing channels. We shall demonstrate in Sec. 4 that for a wide 
choice of model parameters magnitudes of A' relevant for models of rui, are expected to be 
sensitive to the data. 

Once the LHC is in operation the signal size as well as the abihty to probe smaller BR 
are expected to increase dramatically. The task of reconstructing rrif^ , and delineating the 
allowed/disfavored regions of the parameter space in specific models will be much easier. 
The programme for a comprehensive top squark search will certainly take some time. Yet, 
it is gratifying to note that a systematic, largely model independent strategy for top squark 
search in models of mj, is quite possible in a not too distant future. 
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3. The Limiting Values of BR(ti e^d) for Observable 
Signals at Run-II 

In hadron colliders, top squark pairs are produced via gluon-gluon fusion and quark-antiquark 
annihilation, 

gg, qq htl (10) 

The production cross section depends only on the mass of ii without any dependence on 
the mixing angle in the top squark sector, since it is a pure QCD process [40]. The total 
pair production cross section at the Tevatron for y/s = 2 TcV is ~ 15-0.3 pb which is 
40% larger than the cross section for ^/s =1.8 TeV, for the range of rrif^ ~ 100-200. The 
QCD corrections enhance this cross section by ~30% over most of SUSY parameter space 
accessible at Tevatron [31]. 

We investigate the signal of top squark pair production in the channel e~^e~ plus two or 
more jets, assuming that both the top squark decays via a single RPV coupling Ai3j, 

fi ^ e+ + d ; fi ^ e~ + (11) 

where we have suppressed the generation index of the d type quark since we have not 
employed any specific fiavor tagging. 

The leading SM backgrounds corresponding to the signal with opposite sign di- electron 
( OSDE) plus 2 or more jets are the following: 

a. Drell-Yan process via qq' e+e~. 

b. W boson pair production, qq' WW, where both the W decay leptonically, W — > eVe- 
Note that we also consider W decays to r leptons which may decay to electrons. 

c. qq' WZ, where W decays hadronically and Z decays leptonically. 

d. Z boson pair production also leads to the same final state : qq' ZZ {qq') (ee) . 

e. Top quark pair production, qq, gg ^ tt , where both the top quarks decay semi- 
leptonically via W, t ^ beve- 

f. Single top quark production, qq' tb, where one lepton comes from top quark and 
the other comes from 6-quark decay. 

In processes a) - d) additional jets come from initial/final state QCD radiation(I/FSR). 
We have analyzed the signal and background processes using the PYTHIA (v6.206) event 
generator [21]. We generate signal events in the di-electron + jets channel forcing ti to 
decay, ti — > e + g with 100% branching ratio switching off all other allowed decay modes of 
fi in PYTHIA. 

In our calculation we set the rcnormalization and factorization scale to — s and 
CTEQ3L [41] for the parton distribution functions. For the jet reconstruction we use the 
routine PYCELL in PYTHIA [21]. We selected events in the hadronic and electromagnetic 
calorimeter cells in pseudorapidity (rj) and azimuthal angle(0) of size Arj x A0 = 0.1 x 0.1. 
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Cells with Et > 1 GeV are taken as initial seeds to form calorimetric towers. Jets are 
reconstructed with cone radius 0.5 and only those are accepted which has transverse energy 
Et > S GeV and are smeared by 0.5xy/E^. We selected events applying the following set 
of cuts. 

1. Leptons, required to be of opposite charges and of same flavor, are selected with pfp >10 
GeV and \rie\ <2.5 

2. Number of jets arc required to be rij > 2, where, jets arc selected if Ei^ > 15 GeV, 
\r]j\ <3.0. Isolation between any two jets are ensured by demanding AR{j,j) > 0.5, 
where AR = ^/A(f)'^ + Arj'^. 

3. Electrons and jets are assumed to be isolated, if AR{i,j) > 0.5. 

4. Events with di-electron invariant mass between 80 GeV< Mu < 100 GeV and < 
10 GeV are not accepted. 

5. Azimuthal angle between two leptons are required to be (/){££) < 150 degree. 

6. Events are vetoed out for > 25 GeV. 

7. The total visible energy of any event are required to be, St > 350 GeV, where St — 
Hj, + i?^; h!^^^ = scalar sum of transverse energy of all leptons (jets). 

8. We constructed two lepton-jet invariant masses considering all possible combinations of 
the final state particles. Finally we select only that combination in which the difference 
between two is minimum provided \m{iiji) — m{i2j2)\ < 20 GeV. 

Cuts 1-3 are basically event selection cuts. Cut 4 is used to suppress the backgrounds 
where lepton pair is coming due to Z ^ £i decay, where as cut 5 is applied aiming to 
suppress the background due to the Drell-Yan process (a), where leptons are mostly back to 
back in the azimuthal plane. Note that the signal is almost free from any missing momentum 
^. Therefore, using cut 6 we vetoed out those events which involve large amount of missing 
momentum. The background from WW and tt suffer heavily because of this cut. Finally, the 
cut on total visible energy, St > 350 GeV significantly reduces all backgrounds particularly 
DY, to a negligible level, without costing too much in the signal cross section, except for low 
values of rrif^ . 

In Table 1, we summarize our results for all backgrounds. The second column contains 
the raw production cross sections corresponding to each process. In the 3rd column, we 
present the number of events (A^i-e) surviving after cuts 1-6. The effect of cut 7 is shown 
separately in the 4th column(A^7). In the 5th column we present the acceptance efficiencies 
for each of the respective processes. We notice that the jet selection cuts are very effective 
in eliminating the backgrounds due to gauge boson pair productions, as jets are not very 
hard in these processes. In the 1^1^ case, jets mainly arise due to ISR and are very soft. We 

^Some amount of missing momentum may arise if jets or leptons go undetected due to loss in the beam 
pipe, for very low energies which are below the detection threshold or due to energy momentum misman- 
agement. In any case, the missing momentum is not so hard. 
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notice that in this case the selection efficiency turns out to be at the level of ~ 10"*^ due to 
the jet selection cuts. On the other hand, in ZZ and WZ case, the lepton pair comes from 
Z decay, Z — > ee, where as the accompanying gauge boson decays hadronically. Therefore, 
although the jet selection cuts are less stringent, cut 4 and 5 are very effective. Finally, the 
cut on total visible transverse energy drastically reduce all background processes bringing 
them to a negligible level. This table clearly shows that our signal cross section is almost 
background free. The last criterion 8 is used to reconstruct top squark masses and to reveal 
the lepton number violating nature of the underlying interactions. 



Process 


Cross section 


iVl-6 


Nj 


Efficiency 


Events for 




(Pb) 








£=2 fb-^ 


WW 


8.06 


21 





0.0 


0.0 


WZ 


2.34 


487 


5 


1.0x10-5 


0.047 


ZZ 


1.08 


1838 


46 


9.2x10-5 


0.199 


tt 


4.38 


743 


30 


3x10-5 


0.263 


tq 


2.39 


6 





0.0 


0.0 


DY 


2.959x10^ 


91 


1 


5x10-''^ 


2.959 



Table 1: Results of the PYTHIA [21] simulation for all background processes. 

In Table 2 which is of the same structure as Table 1 except for the last two columns, 
we show the signal cross sections for various top squark masses. It is to be noted that the 
cut St >350 GeV costs signal cross section heavily ( by about factors of ~ 10-60) for lower 
values of mf^(~ 80-100 GeV) as leptons and jets are relatively soft where as for higher mf^ 
this cut does not affect the signal cross section significantly. The signal efficiencies vary from 
~ 2-30% for the range of m^- 100-240 GeV. In column 7 we present the significance of the 
signal for e^r =1. The last column presents the minimum value of e^j. that can be measured 
at 5(7 for an integrated luminosity of 2 fb~^. 

In Fig.l, wc show these minimum BRs as a function m^^. The upper region of the solid 
line can be explored by Run-II for C =2 fb~^ luminosity. In the same plane we show the 
region ( above the dashed line) which is already excluded at 95% C.L by Tcvatron data [20]. 
Comparing the two regions we find that the improvement in sensitivity is by ~ factor of 2-3 
for 80 ^ mi^ ^ 160. For higher top squark masses it is still quite significant. As discussed 
in Sec. 2 this is the first step for obtaining a model independent search in the framework of 
RPV MSSM. 

The actual limiting BR may be even smaller as can be seen, a) by replacing the cross- 
sections from PYTHIA ( the second column of the table 2) by the corresponding NLO cross- 
sections of [31] which are typically larger by 30%, b) if accumulated luminosity significantly 
larger than 'iflT^ is considered. Our results arc therefore very conservative. More opti- 
mistic results can be easily obtained by dividing the limiting BR in Table 1 by {cnlo/'^p)^^'^ 
{C_A/2fb~^y/^, where apfio is the cross-section in [31], ap is the PYTHIA cross-section and 
jC^ is the actual luminosity. 

We have not tagged the fiavor of any jet in the final state. We have checked that for 
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140 


1.34 


111196 


28015 


0.0933 


250.0 


134.26 


19.3 


160 


0.617 


111424 


46985 


0.1566 


193.2 


103.76 


21.9 


180 


0.304 


111588 


68625 


0.2287 


139.0 


74.66 


25.8 


200 


0.158 


108769 


83943 


0.2798 


88.41 


47.47 


32.4 


220 


0.084 


106124 


92213 


0.3073 


51.62 


27.72 


42.5 


240 


0.046 


103818 


96204 


0.3206 


29.49 


15.83 


56.2 


260 


0.026 


93410 


89399 


0.2979 


15.49 


8.32 


77.5 



Table 2: Results of the PYTHIA [21] simulation for the signal di-electron plus 2 or more jets 
due to the top squark pair production at Tevatron 

771^^=120(180) the overall efficiency, in Table 2 is reduced to 0.021(0.097) ( including a b- 
tagging efficiency of 50%), if dj is identified with a b-quark. This suppression, however, 
will be adequately compensated by strong reduction in the backgrounds. For example the 
Drell-Yan background will now be non- vanishing mainly due to misidentification of light 
quark and gluon jets as b-jets, the probability of which is extremely small. Assuming the 
signal to be essentially background free and requiring 10 events as the criterion for discovery, 
the limiting BR is found to be 27.3%(41.2%) for mj^ = 120(180). Due to the uncertainties in 
cross-section and £_4 ( see above) these limiting BR may be smaller. We therefore feel that 
the numbers in Table 2 are fairly representative for all d-type flavors. 

As we mentioned in the previous section, the limiting values of ebr can also lead to 
constraints in the MSSM parameter space in RPV models of niy. We discuss them details 
in the next section. In these models the couplings A^33, 7=1-3, are the most relevant ones in 
most scenarios. Considering di-leptons of the same flavor the BR in Fig.l may be interpreted 
as BR (ti eh) or BR(fi pib). 

It is expected that the invariant mass of the lepton and jet should show up a peak at m^^ . 
However, a combinatorial problem arises when the decay of a pair of ti is considered. The 
last kinematic selection 8 is used to reconstruct the top squark mass. The correct lepton-jct 
combination can be separated out by demanding the difference between any two lepton-jet 
invariant mass(777^j) be the minimum. In Fig. 2, we show the lepton-jet invariant mass distri- 
bution normalized for jC—2 fb^^ and with ebr=^ and for three m^^ masses, 100 GeV, 120 GeV 
and 140 GeV which are presented by solid, dashed and short dashed lines respectively. We 
have not shown the corresponding distributions for any of the backgrounds since after im- 
posing all cuts those are turn out to be negligible (see Table 1). As expected, visible peaks 
at each 777^^ is present which are not expected in any of the backgrounds. Therefore, in this 
channel, the mass of ti can be determined with reasonable accuracy. More importantly the 
successful reconstruction of the top squark mass unambiguously imply the lepton number 
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Figure 1: The discovery region above the sohd hne for C =2 fb ^ at Run-II. The region 
above the dashed hne is excluded by Tevatron Run-I data [20]. 

violating nature of the interaction underlying t\ decays. The actual signal size may be con- 
siderably larger due to reasons discussed above. Thus the possibility that the reach will 
extend to higher mi^ or smaller e^^ is therefore quite open. 

4. Stop decay branching ratios and the Hmits on A' in 
models of rriy 

As mentioned in the introduction when the t\ is the NLSP in RPV models of m^^ three decay 
channels are allowed, which may naturally compete with each other in various regions of the 
MSSM parameter space. They are the loop induced flavor changing decay mode, Eq.(5) [22] 
the 4-body decay into states with nearly massless fermions, the bottom quark and the LSP, 
Eq.(6) [27, 35] and the RPV decay mode Eq.(l). In this section we discuss the systematics of 
MSSM parameter space which enable us to identify the regions where different decay modes 
dominate. 

If the sleptons are lighter than ti, then the 3-body decay mode, Eq.(3), involving sleptons 
open up. The competition between Ix hlv and t\ — > hlv and the RPV mode has been 
discussed in ref . [19] . Here we shall also identify regions of the parameter space where the 
decay modes given by Eq.(l), (3) and (5) compete with each other. In the process we also 
demarcate the difficult regions for Run-II in the context of top squark searches, where all 
decay modes may have relatively low rates. 

It has been mentioned earlier that the couplings A^gg arc the most important ones in 
models of mv In [16] the upperbounds on these couplings were obtained from neutrino 
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Figure 2: The lepton-jet invariant mass(m^j) distribution for £ =2 fb ^, for three top squark 
masses 100 GeV(sohd hne), 120 GeV(dashed hne) and 140 GeV(dotted hne). 

data in a variety of scenarii. In all the cases the bounds were found to be approximately of 
the same order of magnitude (~ 10"^ — 10""^). Our analysis based on the limiting BR of the 
last section can estimate the constraints on \'^^- or X^^^- that may arise from Run-II data. 
As before we shall assume only one of these couplings to be dominating and shall henceforth 
drop the index of A'. 

All masses and mass parameters in this paper are in GeV. For our analysis we fix the 
parameters: (i) The CP-odd neutral higgs mass Myi=300, which is required to calculate the 
e parameter (see Eq.(12)), (ii) the trilinear coupling in the sbottom sector ^1^=300 and (iii) 
the trihnear coupling in the stau sector Ai-=200. The variation of the BR with respect to 
the other parameters will be explicitly discussed as and when required. 
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4.1 Competition between the Loop induced and RPV decays 



As is well-known the loop decay width is controlled by the parameter e which denotes the 
amount of iL,R-CL mixing [22] and enters in the decay width, 



where / is the composition of neutralino mixing. The detailed expressions for e and the 
function / can be found in Ref. [22, 27]. Neglecting the lepton masses the decay width of the 
channel ti I'^d ( l—e or jj) is 



where A' is the dominant RPV coupling and 9i is the mixing angle in the top squark sector. 

As long as the 2-body and 3-body RFC decay modes Eq.(2) and Eq.(3) do not open 
up, i.e. if the top squark is the NLSP, the above two modes compete with each other. In 
principle the 4-body decay mode could also enter into the competition. However, in order to 
study the simplest example of competing modes the latter has been suppressed by considering 
relatively large values of tan (5. However, the competition among all these three decay modes 
will be considered later. 

If A' is close to its current experimental bound from indirect searches [17], then the RPV 
decay dominates over the loop decay for the entire region of the parameter space unless cos 9i 
is fine tuned to be very small. The competition between the two modes becomes generic 
when A' is ~ 10~^ — 10~^, which is interesting from the point of view of RPV models of 
neutrino mass [15]. The estimate A' ~ 10~^, as mentioned in the introduction, is based on 
the assumption that the SUSY breaking scale {Msusy) ~ 100 [16]. Somewhat larger values 
of Msusy push this estimate upwards. On the other hand value of A' somewhat smaller than 
~ 10~^ may be relevant if the absolute values of the neutrino masses, which are not known 
at the moment, are much smaller than the typical choice ~ 1 eV. 

As M2 - SU(2) gaugino mass parameter(gaugino mass unification is assumed), // - the 
higgsino mass parameter and tan [3 - the ratio of two vacuum expectation values of higgs 
sector, completely describe the neutralino and the chargino sector. We have chosen these 
parameters such that the m-± is around 200, which more or less fixes the limit of the top 
squark mass up to which the competition among various decay channel can occur. Otherwise 
the 2 body decay mode of ti, Eq.( 2) will be open up and dominate over all other decay 
modes. The common slepton mass is taken to be heavier than mi^ to avoid the 3-body decay 
channel. 

In Fig. 3 the competition between these two decay mode has been illustrated for various 
values of mi^. The other MSSM parameters which are involved in this calculation are: 
M2=250, ji — -1-250, tan /3=40, the common scalar squark mass m^=300, common slepton 
mass m^-=235, cos^j=0.7 and A'=0.001. 

As the top squark mass is increased, the e parameter as well as the phase space factor 

( "s" \ ^ 

1 — in Eq.( 12) increase, but the former rises more sharply. So, although both the 




(12) 




— — A'^mf, cos^r 



(13) 
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widths in Eq.( 12) and Eq.( 13) have a common hnear dependence on nii^, the loop decay 
BR dominates over that of the RPV decay above a certain m^^. This happens for almost all 
choices of the other parameters, unless they are fine tuned to make e very small. 

For smaller value of cos 61, both the loop and RPV decay widths decrease, the former 
through the e term and the latter through the direct dependence on cos^^ respectively. 
The competition between the two BR still occur albeit for higher top squark masses. The 
competition ceases to exist only if cos 6^ is fine tuned to make the e parameter negligible. 




110 120 130 140 150 160 170 180 190 200 

mf, (GeV) 

Figure 3: The RPV and loop decay BRs as functions of m^_^. The Other MSSM parameters are 
M2=250 , n = +250 , tan/?=40, mg=300 , m^-=235 , cos%=0.7 and A'=0.001. 

The RPV decay width depends on the product of cos 9f and A'. Keeping this product ( i.e. 
the width in Eq.(13)) fixed, if we increase A', the loop width will decrease as a consequence 
of lowering cos%. So, the competition will take place for higher top squark masses only. 
However, above a certain A' the loop decay fails to compete for the entire range of m^-^ 
corresponding to a top squark NLSP. On the other hand for smaller A' the RPV decay width 
is scaled down in a straight forward way. Now the competition occurs over a larger range of 
rrif^ and for smaller values of cos^f and /or tan/3. 

If tan f3 is lowered, for fixed fi and M2, the chargino mass is lowered by a small amount so 
that the threshold for the 2-body decay is slightly lowered. More importantly, the e parameter 
decreases dramatically below a certain tan/3. Here the RPV decay overwhelms the loop 
decay. However, precisely for such low values of tan (3 the 4-body decay become important 
if, in addition, the chargino is of low virtuality (wf^ ~ m ~± ; see the next subsection) 

Even if A' is as low as ~ 10~^ the competition between the two modes still exists for 
smaller values of tan (3 which lowers e and. hence, the loop decay width. The minimum value 
BR, as shown in Fig.l can be traded to find the limiting value of A' considering Eq.( 12) 
and 13. In Fig. 4 the two curves represent limiting values of A' for observable signal for 
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Figure 4: The minimum value of A' at 5a for observable signal for cos 0^=0.7 (solid) and 0.02(dashed 
curve), the others parameters are the same as in Fig. 3 

two values of cos The regions above the curves corresponds to observable BR as given in 
Fig.l. The other SUSY parameters chosen are as in Fig. 3. In this figure and the similar 
ones presented subsequently the horizontal arrow represents the upper bound on X[^i prior 
to the neutrino data [17]. The bounds on and A23J are even weaker. Only the bound on 
A^33 is ~ 10""^. Hence significant improvement in the existing limits on many RPV couplings 
is expected. For larger cosO^, the RPV width increase significantly. As a result the BR 
constraint is satisfied for lower A'. The sharp rise in the curve for rrif^ ^200 is a consequence 
of the opening up of the 2-body channel ii bxf- It is interesting to note that for large 
cosOj: the data will be sensitive to the values of A' relevant for neutrino masses until the 
2-body decay channel opens up. 

4.2 Competition between the 4-body and RPV decay 

The dependence of the 4-body decay rate on supersymmetric parameters has been discussed 
in great detail in Ref. [27, 35]. The competition between 4-body decay modes with the loop 
induced fiavor changing decay mode Eq.(5) has been discussed both in MSSM and mSUGRA 
models in Ref. [35]. 

In general in order to identify the parameter space relevant for the competition between 
the RPV decay channel and the 4-body decay channel, one has to take the ti to be almost 
right handed (i.e., cos 9^ small) and A' ~ 10~^ or 10""^. For small value of tan /3, the loop decay 
amplitude becomes neghgible without requiring any fine tuning. In Fig. 5 we demonstrate 
this competition for A' = 10"^ as a function of m^^. The choice of the MSSM parameters are 
explicitly mentioned in the figure caption. 
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Figure 5: The RPV and 4-body decay BRs are shown. The Other MSSM parameters are M2=250 
, II = +250 , tan/3=6, mg=300 , m^-=210 , cos6'£=0.1 and A'=0.0001. 




120 140 160 180 200 220 240 260 

m,- (GeV) 

Figure 6: The minimum value of A' at 5(T for observable signal for cos 0^=0. 7(solid), O.l(dashed) 
and 0.02 ( dotted curve). The other parameters are the same as in Fig. 5 

As expected the 4-body decay channel opens up for relatively low (m^^ - m-± ) so that the 
chargino in the 4-body decay process has a small virtuality. As before, in Fig. 6 we show the 
potential region(upper side) of A' which can be probed at Run-II experiments. This region 
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Figure 7: The RPV, loop and 4-body decay BRs are shown. The other MSSM parameters are 
M2=250 , II = +250 , tan/?=10, mg=300 , m^-=210 , cos6'£=0.9 and A'=0.0001. 

corresponds to the BR constraint (Fig. 1) are shown for three values of cos % and other SUSY 
parameters chosen are as in Fig. 5. The conventions are same as in Fig. 4. 

4.3 Competition between the Loop, 4-body and RPV decay 

In order to illustrate the possibility of competition between the above three channels, we 
shall keep in mind that the e parameter must not be as small as in the previous section. 
The competition is demonstrated in Fig.7 with the choice of SUSY parameters mentioned 
explicitly in the figure caption. 

The total 4-body BRs is significant (~ 10%) for the range of top squark mass, m^± - 20 
~ m^^ S; m -± . In this top squark mass range the chargino is of very small virtuality. Also 
we choose i.e. the common slepton mass, such that even after the mixing in the third 
generation, the lighter tau slepton mass m^^ is above the chargino mass. So the slepton 
mediated 4-body process also has a low virtuality. But the 3-body decay mode Eq.( 3) is 
still kinematically forbidden. 

If the signal is seen in all three channels then one has to identify the region of the 
parameter space where the corresponding BRs are above the observable limit. Similarly in 
order to exclude a particular m^^ comprehensively it is essential to establish that at least 
one of the competing modes would be observable over the entire parameter space. In order 
to do a complete job one needs the minimum observable BRs at Run-II for each of the 
allowed modes. Unfortunately at the moment we have numerical estimates for the RPV 
mode (ti — > e'^dj) only ( Fig.l). In the following we shall delineate the regions of the 
parameter space where i) the RPV decay is observable at Run-II or ii) either one of the two 
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Figure 8: Different regions dominated by a particular decay channel of top squark are shown for 
m£j=180 . See the text for conventions for demarcating regions. Except for cos 6*^=0.3, all the other 
parameters are the same as in Fig. 7. 
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Figure 9: Similar information as is in Fig. 8 using different variables. 



competing modes have sizable BR. 

The relevant information will be presented in the form of scatter plots obtained by varying 
two important parameters randomly keeping the others fixed. The scatter plots also illustrate 
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m,- (GeV) 

Figure 10: Similar information as in Fig. 6, for cos 0.9 (solid), O.l(dashed) and 0.02(dotted 
curve). Other parameters are the same as in Fig. 7 

the competition among the decay modes in specific regions of the parameter space. 

In Fig. 8, fixing mf^=180, cos^j=0.3, and tan/3 are varied randomly setting the other 
parameters as in Fig. 7. The fixing of m^^, which tacitly assumes that rrif^ can be recon- 
structed, make the analysis simpler. The systematic of the parameter space is clear from 
Fig. 8. In this figure, the regions marked by the circle is the one where the RPV mode is 
above the observable limit, i.e., BR( ii ^ + d) ^ 26%(see Fig.l). Although the width 
of this mode does not depend upon // or tan/3 directly, however, its BR is quite sensitive 
to these parameters. The regions marked by '+ ' correspond to BR( ti cxi)~ 75% and 
that marked by the black diamonds correspond to BR {ti ~ 30% with the RPV 

BR less than the observable limit. Note that for low values of tan/5, the RPV BR is much 
larger than 30%. The region labeled by 'A', rrif^ > m^± +mh, the 2-body decay mode Eq.(2) 
opens up and overwhelms all other decay channels. Finally, the region marked by 'B', where 
/X and tan/3 are large, the lighter r slepton mass eigenstate (fi) becomes rather light and 
the 3-body decay mode involving a fi in the final state strongly dominates. In the dotted 
region fi is lighter than the Xi or has unphysical mass. Although m^^ ~ m^o is allowed in 
RPV MSSM in general, but we are not interested in top squark signals in this scenario. 

If the top squark signal is seen in one or more channels then one can broadly identify 
the relevant region of the parameter space. For example if all the three modes are seen then 
the white region ('C') or regions in its neighborhood may be of interest. For more precise 
conclusion one needs to know the limiting BR of all the modes quantitatively . One can hope 
that Tevatron Run-II and/or LHC will gradually supply the relevant information. However, 
the same region may eventually turn out to be the difficult one to exclude at Run-II, if 
no signal is seen, since in parts of this region all the BRs may turn out to be below the 
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Figure 11: The RPV, loop and 3-body decay BRs as funetions of m^^. The other MSSM parameters 
are M2=250 , = +250 , tan/?=40, mg=300 , m^-=175 , cos6'j=0.5 and A'=0.001. 

observable limit. 

In Fig. 9 the scatter plot is in tan /3 — cos 9i plane with //=250 and the other parameters 
as in Fig. 7. The convention for demarcating the regions arc the same as in Fig. 8. Again the 
white regions could be the difficult ones from the point of view of comprehensive t\ search 
at Run- 11. 

In Fig. 10 the three curves represent limiting values of A' for observable signal for three 
values of cos 9i. The regions above the curves corresponds to observable BR as given in Fig.l. 
The other SUSY parameters chosen are as in Fig. 7. It is interesting to note that for large 
cos 9i the data will be sensitive to values of A' relevant for neutrino masses until the 2-body 
decay channel opens up. It is to be noted that for relatively small mi, the bound is fairly 
insensitive to cosOi for the range 0.1- cos^- 0.9. As the threshold of the 4-body decay open 
up for larger mi, , the constrain on A' gets weaker as expected. 

4.4 Competition between the 3-body, Loop and RPV decay 

The competition between RPV decay mode, the loop decay, and all RFC 3-body channels 
has been studied in [39]. In this section we consider a scenario where top squark is not 
the NLSP and the first two RFC decay modes of Eq.( 3) are open. Wc have then studied 
the competition among these two modes, the loop decay and the RPV decay taking into 
account the limiting BR of the last mode obtained in Sec. 3. As the 3-body decay mode is 
kinematically allowed for light sleptons only, the slepton mass should be chosen with care so 
that it is consistent with the experimental lower limit. 

With the choice of the SUSY parameters as in Fig. 11, the 3-body decays, if kinematically 
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Figure 12: Similar information as in Fig.8, but for mj^=160 . The other parameters are the same 
as in Fig. 11. 
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Figure 13: Similar information as in Fig.l2 using different variables. 

allowed, have BR ~ 10% almost for the entire range of top squark masses. For this set of 
parameters, the m^^, mj>, and m^-^ (/ = e or yu) are 124 , 156 and 175 respectively. 

Interestingly we have found that when the chargino is in the mixed region with a relatively 
large mass, i.e., when the 3-body decay width is somewhat reduced both due to a mixing 
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Figure 14: Similar information as in Fig.6 for cos 6*^=0. 5(solid), O.l(dashed) and 0.02(dotted curve) 
and the others parameters are the same as in Fig. 11 

angle factor and propagator suppression, there may be a competition among the loop, 3-body 
decay and RPV decay modes for A'=0.001. Here the loop decay width is significant thanks 
to relatively large cos^f=0.5 and large tan/9. This is demonstrated in Fig. 11. It follows 
from Fig. 11 that in the neighborhood of mi^ = 150 all three modes coexist with nearly equal 
BR. However from the limiting BR plot( see Fig.l ) we find that in this region the signal 
is observable if BR(ii — > e+ + d ) ^20%. Hence the loop decay channel may not be very 
important as the discovery channel. If the BR(fi — > e+ + d ) is below the observable limit, 
the 3-body mode will be the main discovery channel. In Fig. 11 for relatively low m^^ the 
3-body mode with fi in the final state opens up. For the higher m^^ the modes with u and 
other sleptons in the final state are allowed. 

In Fig. 12 the competition among the decay modes is illustrated in /i — tan/? plane, for 
mf^ =160 . For this mj^ only the 3-body mode with the fi in the final state is relevant. 
Depending on the values of tan /3, and their product in different regions of the parameter 
space the competing decay modes are kinematically allowed. In Fig. 12 the dotted circles 
delineate the parameter space where the RPV decay is observable. The regions characterized 
by relatively large (j, and tan (3 corresponds to the light fi scenario. This part of the parameter 
space dominated by the 3-body decays ( BR ^ 70% is marked with the black diamonds. 
The dotted region is theoretically disfavored as explained in the context of Fig. 8. Finally, 
the black circles represent the parameter space with 40% ~ BR(ti cxi)^70%. Only a few 
points appear at large tan/3. In the region marked with 'A' the 2-body decay overwhelms 
the other modes. In Fig. 13 a tan/3 — cos% scatter plot is presented following the same 
convention. 

In Fig. 14 we present the hmiting value of A' for three values of cos^f corresponding to 
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the parameter space of Fig. 11. The first change in the slope occurs at about mf^=13Q due 
to the opening up of the channel ti — > buT-fi. The second change in the neighborhood of 
mf^=180 corresponds to the decay mode ii — > buih { I — e or ji) . In both cases the BR 
of the RPV decay mode arc reduced which have to be compensated by higher values of A'. 
Finally for m^^ ~ 210, the 2-body channel, Eq.(2), becomes the main decay mode. 

5. Conclusion 

It is quite possible that the mass of the lighter top squark is much smaller than the other 
squarks and gluinos due to mixing and RG effects and it is the only strongly interacting 
superparticle within the kinematic reach of Run-II of the Tevatron with large production 
cross-section. If this is the case then the RPV decay ti ifdj driven by the trilinear 
coupling A-3j, where i and j are generation indices, may be the most attractive channel for 
discovering R-parity violation[18, 19, 20]. 

Additional interest in this process stems from the fact that some subset of the above 
couplings, in particular A^gg, may be important ingredients of RPV models of m,^[15, 16]. This 
scenario constraints the magnitudes of these couplings to be generically small (~ 10~^ — 10""^, 
see, e.g., [16]). 

If the couplings are indeed so small the RPC 2-body decay, (Eq.(2)) or 3-body decay 
modes (Eq.(3)), if kinematically allowed, would overwhelm the RPV decay and the LSP 
decay may be the only signature of R-parity violation [25]. This signature, Eq.(4), however, 
may not reveal the lepton number violating nature of the underlying interaction or whether 
the strength of the coupling is indeed in the right ballpark required by models of rriiy. 

The situation is dramatically different if the ti is the NLSP since the allowed RPC decays 
- the loop induced ( Eq.(5)) or the 4-body ( Eq.(6)) channel - are naturally suppressed. If 
the RPV couphng is indeed ~ 10~^ — 10"^ then BR of the three allowed channels may indeed 
be comparable. Thus the simultaneous observation of two or more of these decay may be a 
hallmark of RPV models of rriu. 

In Sec. 3 using event generator PYTHIA[21] we have estimated the minimum value of BR 
of the RPV decay channel ii — > e'^dj for various values of mf^ corresponding to observable 
signals at Run-II experiments. Our results (see Fig.l) show that much smaller BR can be 
probed at Run-II with 2fb^^ of data compared to the bounds obtained from Run-I data 
[20]. These results are approximately valid for down type quarks of all generations and 
also for the channel ti — > ii^dj. In reality the limiting BR may be much smaller than 
our conservative estimates as can be seen by using enhanced NLO cross-section[31], larger 
integrated luminosity or by employing b-tagging, since in many models A^gg are the most 
important couplings, to improve the S/\/B ratio. Our simulations show that m^^ can be 
reconstructed from the decay products with reasonable accuracy, revealing thereby the lepton 
number violating nature of the underlying decay dynamics. 

It is gratifying to note that even our conservative estimates of the limiting BR can be 
translated into interesting upper bounds on the RPV couplings X[^j {i—e or //) for represen- 
tative choices of the MSSM parameters if no signal is seen (Figs. 4, 6, 10,14). Thus the 
existing bounds [17] on several A^g - ( except perhaps A^gg ) can be significantly improved. 
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These results indicate that the Run-II data will indeed be sensitive to magnitudes of these 
couplings even if they are as small as that required by the models of m,^. 

Using our estimate of the limiting BR as a function of mf^ one can demarcate the regions 
of the MSSM parameter space in specific models, where the RPV decay is observable. In Sec. 
4 we have also studied the systcmatics of the MSSM parameter space and have delineated 
the regions where the competing decay modes are numerically significant. One can also 
have some idea of the difficult regions of the parameter space where the BR of none of 
the competing decays clearly dominates. All these information will become more precise 
once full simulations of the competing signals (Eq.(7) and Eq.(8)) estimate the limiting BR 
corresponding to all signals . If no signal is seen then programme for top squark search at 
the LHC may focus on the regions of the parameter space, which were difficult in Run-II 
experiments. 
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